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The Asteroid Robotic Redirect Mission (ARRM) has been the topic of many mission 
design studies since 201 1. 1 The reference ARRM spacecraft uses a powerful solar electric 
propulsion (SEP) system and a bag device to capture a small asteroid from an Earth-like 
orbit and redirect it to a distant retrograde orbit (DRO) around the moon. The ARRM 
Option B spacecraft uses the same propulsion system and multi-Degree of Freedom (DoF) 
manipulators device to retrieve a very large sample (thousands of kilograms) from a 100+ 
meter diameter farther-away Near Earth Asteroid (NEA). This study will demonstrate that 
the ARRM Option B spacecraft design can also be used to return samples from Mars and 
its moons - either by acquiring a large rock from the surface of Phobos or Deimos, and/or 
by rendezvousing with a sample-return spacecraft launched from the surface of Mars. 


I. Introduction 

The Asteroid Robotic Redirect Mission (ARRM) Option A concept, first introduced in 2011, 1 is a mission 
design to capture and return a small near-Earth asteroid (NEA) to cislunar space. The ARRM Option B 
concept is a similar spacecraft, but designed to return a very large sample from a more difficult to reach 
NEA. 2 In this work we show that the spacecraft designed for ARRM Option B is also well-suited to sample 
return from Mars and its moons. This work presents low-thrust interplanetary trajectories from cislunar 
space to Mars and back, including descent from the Martian sphere of influence (SOI) to the desired orbit 
altitude and ascent to the SOI after the sample retrieval is complete. Missions are presented to Phobos, 
Deimos, and low-Mars orbit (LMO), the latter from which the spacecraft could retrieve a sample capsule 
launched from the surface of Mars. 

Boulder returns from Deimos or Phobos, or a Mars sample return are of notable value. The two moons 
are considered potential bases for manned teleoperations outposts, and serve as stepping stones to eventual 
human exploration of Mars (boots on ground). 3,4 A precursor mission substantially minimizes the risk of 
any subsequent human mission to these bodies, enabling an examination of their composition. A robotic 
return of moon material would also help to answer questions as to their origin. Similarly, a Mars sample 
return has been charted as a critical step on the path to placing humans on the planet. It is the highest 
priority flagship mission identified in the 2013-2022 planetary science decadal survey. 5 Asteroid robotic 
redirect vehicle (ARRV) trajectory performance is examined for viability of rendezvousing with a sample in 
a 400 km altitude orbit that is launched from the surface by another system. 

Mission design is carried out using the Evolutionary Mission Trajectory Generator (EMTG), NASA 
Goddard Space Flight Center (GSFC)’s low-thrust trajectory optimization tool. 6,7,8,9,10,11,12 This work 
presents the Mars System Sample Return (MSSR) problem as a global optimization problem in the EMTG 
framework and models the complex relationship between the mass of the sample, the mission time constraints, 
and the trajectory. EMTG is used to design the interplanetary phases of the trajectory and the Edelbaum 
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approximation is used to model the flight time and propellant use from the SOI to Deimos/Phobos/LMO 
and back to the SOL 


II. Mission Modeling 


II. A. Mission Architecture 

The MS SR architecture considered this work is based heavily on the ARRM Option B mission concept. 
ARRM Option B consists of a large solar-electric propulsion (SEP) flight module with an attached mission 
module. The combined vehicle is referred to as the asteroid robotic redirect vehicle (ARRV). The mission 
module is made up of a set of gripper arms, a mechanism to anchor to the target, a structure on which to 
mount them, and avionics. In ARRM Option B, the grippers are used to acquire a boulder from the surface 
of a large near-Earth object (NEO) and hold that boulder all the way back to Earth. The MSSR paradigm is 
identical except that the grippers are used either to acquire a boulder from the surface of Phobos or Deimos, 
or to grapple a vehicle launched from the surface of Mars. 

The MSSR spacecraft is powered by a pair of large solar arrays that supplies 51 kW of power at beginning 
of life (BOL) at 1 AU. The performance of the arrays as a function of distance from the sun is modeled 
using EMTG’s power model 10 with a set of coefficients for triple-j unction gallium arsenide (TJGA) cells and 
the array is assumed to lose 1% of its performance for each year of flight. The spacecraft bus consumes a 
constant 0.5 kW and a 15% margin is applied on input power to the propulsion system. 

There are two versions of the MSSR propulsion system, both based on a high-power Hall thruster. The 
first version of the propulsion system is modeled as a single “super-thruster” running at 2000 s specific 
impulse ( I sp ) and 55% system efficiency that accepts up to 40 kW of input power. Thrust is computed 
using EMTG’s propulsion model 10 as a function of available power. The second version is also modeled as 
a “super-thruster” but runs at 3000 s I sp with 60% system efficiency. The 2000 s version is slightly heavier 
than the 3000 s version and adds 150 kg to the spacecraft. Both propulsion systems are composed of the 
same hardware, but the 2000 s system is run at 400 V while the 3000 s system is run at 800 V. There are 
two additional thrusters in the 2000 s system. Both propulsion systems are run at a 90% duty cycle. An 
11% margin is carried on the xenon propellant. 

There are two launch paradigms for MSSR. The first is to launch the spacecraft on the Space Launch 
System (SLS) Block la launch vehicle to a C3 of up to 27 km 2 / s 2 . The SLS can deliver up to 16072 kg to 
this energy based on current notional performance curve. A 30 day coast period is enforced after launch for 
spacecraft systems checkout. The alternative launch paradigm is to launch to trans-lunar injection (TLI) 
on a Falcon Heavy (FH) and then perform a double lunar flyby sequence. FH is assumed to be capable of 
launching 13200 kg to a C3 of -2.0 km 2 / s 2 which enables the lunar flyby sequence. The double flyby delivers 
the spacecraft to a C3 of 2.0 km 2 / s 2 and takes 60 days, which conveniently is enough time to perform systems 
checkout so no coast period is required after the spacecraft escapes from cislunar space at the second flyby. 13 
While the Earth departure portion of the trajectory is critical, it is not the focus of this study and the results 
shown allow for a departure C3 values that are less than or equal to 2.0 km 2 / s 2 . The practicality of this 
assumption has been verified by reproducing the trajectories generated in EMTG in end-to-end simulations 
in the high-fidelity modeling of the Copernicus software. 14 

Upon arrival at Mars, the spacecraft descends to the altitude of Phobos, Deimos, or a 400 km altitude 
LMO depending on the desired mission, spends between 200 and 600 days, and then ascends to the boundary 
of the Mars SOI and then returns to Earth. A constraint is applied to Earth return such that the incoming 
C3 is less than or equal to 2.0 km 2 / s 2 , which allows for a lunar flyby and then entry into a distant retrograde 
orbit (DRO). The Av necessary to insert into the DRO from the arrival trajectory is approximated at 100 
m/s. Finally, the dry mass of the spacecraft is constrained to be no less than 5170 kg for the 2000 s I sp 
version and no less than 5070 kg for the 3000 s I sp version. These parameters are listed concisely in Table 

I. No other information was supplied to EMTG. 

II. B. The Sims-Flanagan Transcription 

The Sims-Flanagan transcription is a widely used method in which the continuous-thrust trajectory is dis- 
cretized into many small time steps, and the thrust applied during each time step is approximated as a small 
impulse placed at the center of the time step. The trajectory is propagated between control points by solving 
Kepler’s problem. 15 The Sims-Flanagan transcription, when used with an nonlinear programming (NLP) 
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Power Assumptions 

BOL Power at 1 AU 

51 kW 

Array performance model 

TJGA 

Array decay rate 

1% per year 

Spacecraft bus power 

0.5 kW 

Power Margin 

Propulsion Assumptions 

15% 

Thruster model 

Fixed I sp and system efficiency 

Input power bounds 

0-40 kW 

Thruster performance 

2000 s I sp with 55% efficiency 
3000 s I sp with 60% efficiency 

Duty cycle 

90% 

Propellant margin 

11 % 

Launch Assumptions 

Maximum departure hyperbolic excess energy (C3) 

2 km 2 /s 2 (FH) 
27 km 2 /s 2 (SLS) 

Maximum launch mass 

13153 kg (FH with 2000 s I 8p ) 
13169 kg (FH with 3000 s I sp ) 
16072 kg (SLS) 

Declination of Launch Asymptote (DLA) bounds 

+/ — 28.5° 

Post-launch forced coast 

30 days (SLS) 
none (FH) 

Mars stay assumptions 

Desired altitude 

20063 km (Phobos) 
5844 km (Deimos) 
400 (LMO) 

Stay time at desired altitude 

Earth Return Assumptions 

200-600 days 

Maximum arrival C3 

2.0 km 2 / s 2 

Post-mission Av 

100 m/s 

Minimum spacecraft final mass 

5170 kg (2000 s I sp case) 
5020 kg (3000 s I sp case) 


Table 1: ARRM Option B Assumptions 
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solver such as Sparse Nonlinear Optimizer (SNOPT) and a suitable initial guess, is very fast and robust. 
It is considered to be a “medium-fidelity” transcription and is used in existing software packages such as 
Mission Analysis Low-Thrust Optimization (M ALTO), 16 Gravity Assisted Low-thrust Local Optimization 
Program (GALLOP), 17 and Parallel Global Multiobjective Optimizer (PaGMO). 18 

In the classical Sims-Flanagan transcription, the optimizer chooses the three components of an impulsive 
Av vector at the center of each time-step. In order to improve the robustness of the solver, a modified 
transcription known as “up-to-unit vector control” is used in EMTG, where instead of choosing the Av vector 
directly the optimizer instead chooses a control 3-vector in [—1.0, 1.0] that is multiplied by the maximum 
Av that the spacecraft can produce in that time-step. The magnitude of the control vector is then bounded 
in the range [0.0, 1.0], i.e., 


A Vj = UiAv max ,i, ||ui|| < 1.0 


( 1 ) 


where 


a i Dn ava ii a bl e T maiX ( tf to) / Q \ 
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where D is the thruster duty cycle, n ava ii a bie is the number of available thrusters, T max is the maximum 
available thrust from one thruster, to and tf are the beginning and ending times of the time step, m is the mass 
of the spacecraft at the center of the time step, and N is the number of time steps in the phase. This modified 
Sims-Flanagan transcription is used in MALTO, Parallel Global Multiobjective Optimizer (PAGMO), and 
EMTG. 

In each phase of the mission, the trajectory is propagated forward from the first endpoint (i.e. planet) 
and backward from the second endpoint. The trajectory is propagated by solving Kepler’s equation and 
the spacecraft mass is propagated by assuming a constant mass flow rate across the each time-step. The 
specific Kepler propagator algorithm used in EMTG is a Laguerre- Conway method. 19,20 A set of nonlinear 
constraints are applied to ensure continuity in the center of the phase, 
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( 3 ) 


The optimizer also chooses the initial and final velocity vectors for each phase. If a phase begins with 
a launch, the magnitude of the initial velocity vector is used with a launch vehicle model to determine the 
initial mass of the spacecraft as described later in this work. If a phase begins with a planetary flyby, 
two nonlinear constraints are applied to ensure that the flyby is feasible. First, the incoming and outgoing 
velocity vectors with respect to the planet must be equal, 

*4 - = 0 ( 4 ) 

where and are the velocities before and after the flyby, respectively. Second, the spacecraft may not 
fly closer to the planet than some user-specified minimum flyby distance: 
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where 


5 = arccos 



(6) 


Here ii p i ane t is the gravitational parameter of the planet, r p i ane t is the radius of the planet, S is the flyby 
turn angle, and h sa f e is the user-defined minimum altitude. 

Figure 1 is a diagram of a simple low-thrust mission to Jupiter with one Earth flyby using the multiple 
gravity assist with low-thrust (MG ALT) model. The continuity constraints are deliberately left unsatisfied 
in the diagram to illustrate where they must be applied. 

There are four significant advantages to using the Sims-Flanagan transcription. First, the optimal ob- 
jective function value for a Sims-Flanagan trajectory design is usually very close to the optimal cost value 
for a higher-fidelity version of the same trajectory. Second, a low-thrust trajectory generated using the 
Sims-Flanagan transcription makes a very good initial guess for a higher-fidelity trajectory design. Third, 
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Figure 1: An Example Trajectory Using the Sims-Flanagan Transcription 


the Sims-Flanagan transcription is very fast because it does not require numerical integration of differential 
equations. Fourth, the convergence of an NLP solver solving a Sims-Flanagan problem is very robust to poor 
initial guesses, making it ideal for an automated design approach. It is therefore very reasonable to use the 
Sims-Flanagan transcription for preliminary design of missions like MS SR. 

II. C. Spiral Low-Thrust Trajectory Modeling via Edelbaum’s Method 

Low thrust orbit raising and lowering about a planetary body is often referred to as “spiraling” because the 
time-scale of the orbit is short relative to the length of the maneuver and so the spacecraft performs many 
revolutions about the planet before arriving at its desired altitude. When the trajectory is viewed it looks 
like a spiral. Because of the large number of revolutions it is impractical to design such a trajectory using the 
Sims-Flanagan transcription because far too many time-steps would be required and therefore the problem 
would become intractible. Fortunately there is a convenient approximation that makes the low-thrust orbit 
raising and lowering problem much easier. 

Edelbaum 21 showed that when the magnitude of the applied thrust is very low relative to gravity then 
the spacecraft’s orbit may be approximated as being circular at all times. For “up spiral” trajectories, thrust 
is applied along the orbit velocity direction until the orbit reaches the boundary of the planets SOI where the 
orbit becomes hyperbolic with an excess velocity of 0. Thrust is applied opposite the orbit velocity direction 
for “down spiral” trajectories. No plane change is necessary during either type of spiral trajectories because 
the excess velocity of 0 can be reached from any orbit plane of choice. The Av required to transfer from 
a higher orbit to a lower orbit is therefore simply the difference between the circular orbit speeds at the 
starting and ending altitudes, i.e., 


A I /^planet I ^planet /«\ 

Edelbaum — \ / \ {• ) 

V ro V r f 

where ro and rf are the radii of the starting and ending orbit, respectively. In EMTG a spiral arrival is 
determined by solving Equation 7 where ro is equal to the radius of the planet’s SOI and rf is the radius of 
the desired final circular orbit. Conversely, for departure spirals ro is the radius of the initial circular orbit 
and rf is the radius of the planet’s SOL 

The propellant necessary to accomplish the spiral may then be determined by, 
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where rribef ore-spiral is the mass of the spacecraft before the spiral orbit transfer and go is the acceleration 
due to gravity at sea level on Earth, 9.80665 m/s 2 . The time required for the spiral may then be computed 
by determining how much time is required for the thruster to expel the necessary propellant, 


^Edelbaum 


TTl propellant 
171 


(9) 


where m is the mass flow rate of the thruster, which is a function of the available electric power according to 
EMTG’s propulsion and power modeling system. 10 In EMTG, the available power and therefore the thrust, 
I sp , and mass flow rate are held constant across the spiral orbit transfer, fixed to the power available at 
the beginning of the maneuver. Therefore for an arrival spiral, EMTG computes the power at the epoch 
where the spacecraft is at the edge of the planet’s SOI and for a departure spiral EMTG computes the power 
available at the epoch where the spacecraft is in its initial circular parking orbit about the planet. 

The Edelbaum approximation, while not an exact model of the spiral descent and ascent problem, is a 
very good approximation suitable for preliminary design. In this work the Edelbaum equations are used to 
model the motion of the spacecraft from the radius of Mars’s SOI to the orbital altitudes of Phobos and 
Deimos, and to LMO, and back up to the SOI boundary. 


III. Trajectory Optimization 

III. A. The Rendezvous Sample-Return Problem with Fixed Pickup Mass 

The first case considered in this work is zeroth-order approximation of the true Mars sample return problem. 
The spacecraft travels to Mars in the heliocentric frame, disappears, and re-emerges some time later with 
its mass incremented by the difference between the size of sample collected and the mass of the propellant 
necessary to perform all maneuvers in the Mars reference frame, i.e., 


^ Mars— departure — ^ Mars— arrival T 171 sample ^propellant— used— at— Mars (10) 

where m prope ii ant - use d-at-Mars is the mass required to descend to and return from the orbit of interest, i.e. 
that of Deimos, Phobos, or LMO. This formulation of the problem is most appropriate if m sarnp i e is varied 
parametrically and m prope u ant -used-at-Mars is calculated external to the trajectory program. The MSSR 
problem then becomes a straightforward two-journey out-and-back mission to Mars which can be designed 
in any of several existing low-thrust trajectory design programs. 10,16,17,18 

The Jacobian sparsity pattern for the rendezvous, fixed pickup mass MSSR problem with no planetary 
flybys is shown in Figure 2. The diagonal blue lines represent the dependence of the throttle constraint 
in each time-step on the control variables in that time- step. The horizontal solid blue bars represent the 
dependence of the phase match-point constraints on all of the control variables in each phase and on time. 
The sparse set of dots at the bottom represent the dependence of the return C3, return date, total flight time, 
and dry mass constraints on the launch date, the initial mass of the spacecraft, the flight times for each phase, 
and the stay time at Mars. The solitary blue dot in the center of the bottom row represents the dependence 
of the objective function, i.e. maximum return mass, on the mass at the beginning of the return journey. 
The pattern in Figure 2 does not take into account any dependencies introduced by the spiral descent and 
ascent because they are computed outside of the main trajectory optimization. Note that, if the off-line 
computation of the Edelbaum equations is omitted, this formulation also works for retrieving a sample from 
an asteroid as in the original ARRM Option B concept. This sparsity pattern is useful for comparison with 
those of other formulations of the MSSR problem later in this work. While the rendezvous, fixed pickup 
mass MSSR problem is convenient to formulate in most existing tools, it comes with a significant drawback. 
A series of nested parametric runs must be performed to explore the trade-space. First, for a variety of 
launch dates, outbound trajectories must be designed from Earth to Mars. An offline batch computation of 
the Edelbaum equations must then be conducted to find the time and propellant mass necessary to descend 
from the boundary of Mars’s SOI to the target orbit. Pickup mass is then varied parametrically and, for each 
pickup mass, spacecraft arrival scenario at Mars, and desired stay time at Mars, a return trajectory must 
be computed. Let Ned be the number of launch dates in the parametric sweep, and let N t _ out b oun d be the 
number of outbound flight time discretizations, Npm be the number of pickup mass discretizations, N stay 
be the number of stay time discretizations, and N return -date be the number of return date discretizations. 
Then, since the time and propellant mass required for the spiral descent and ascent are a function of the 
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Figure 2: Sparsity Pattern for the Rendezvous Sample- Return Problem with Fixed Pickup Mass 


other variables, the number of trajectories to be computed in the parametric sweep is: 


N total A LD-^t— outbound^ PM A stay -^return— date (ii) 

The expression in Equation 11 can quickly grow to an impractically large number of trajectories. It is 
therefore desirable to find a more efficient method to perform the search. 

III.B. The Sample- Return Problem with Fixed Pickup Mass and Spiral Modeling at Mars 

A better, but still not fully modeled, version of the MSSR problem can be constructed by including the 
computation of the Edelbaum (or other low-thrust spiral approximation) equations in the trajectory tool 
itself as is done in EMTG. 10 This approach removes the necessity to perform any off-line computations and 
therefore allows the outbound and return components of the mission to be optimized in a single run, cutting 
down on the total number of runs of the trajectory optimization tool. However, in this formulation pickup 
mass is still varied parametrically and so the number of runs is still large: 


Ntotal = ^LD^PM^return —date (12) 

The Jacobian sparsity pattern for the fixed pickup mass with spiral modeling formulation of the MSSR 
problem with no planetary flybys is shown in Figure 3. The entries which are common between the fixed 
mass, rendezvous problem and the fixed mass, spiral problem are shown in blue. The new entries which are 
introduced by including the spiral descent and ascent in the trajectory optimization are shown in red. The 
single red dot in the bottom left represents the new dependency of the flight time constraint on the launch 
date. This dependency occurs because the launch date influences the arrival date at Mars, which in turn 
influences the distance between Mars and the Sun at arrival, which in turn influences the available power 
for the spiral descent to the desired orbit and therefore the available thrust and I sp . This dependency on 
solar distance influences the amount of time necessary for the spiral descent and ascent and therefore the 
total mission time. Similarly the line of three red dots next to the solitary red dot represent the influence 
of the outbound journey flight time on the Earth return date, the total flight time, and the spacecraft dry 
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mass constraint. Like the launch date, the outbound flight time influences the date of arrival at Mars and 
therefore the available power and thus the flight time and propellant mass necessary to perform the spiral 
descent and ascent. These variables in turn influence the return date and total flight time as expected but 
also the total amount of propellant consumed that affects the dry mass constraint. These new dependencies 
were implicit in the fixed mass, rendezvous version of the MSSR problem but did not need to be considered 
by the optimizer. Switching to the fixed mass, spiral modeling version of the problem allows the total number 
of optimizer runs to be reduced but require that the optimizer properly process the relationship between 
time and mass in the spiral descent and ascent modeling. 



Figure 3: Sparsity Pattern for the Sample- Return Problem with Fixed Pickup Mass and Spiral Modeling 


III.C. The Sample-Return Problem with Variable Pickup Mass and Spiral Modeling at Mars 

While the fixed mass, spiral modeling version of the MSSR problem is simpler to use and can explore the 
trade space in fewer runs of the trajectory optimization program, it is still not ideal because the pickup mass 
must still be discretized. It is desirable to find a formulation of the MSSR problem that does not require 
any discretization at all except, if desired, on launch date and return date. This approach is implemented by 
allowing the optimizer to directly vary the pickup mass as is done in EMTG. 10 There is therefore no need 
for a parametric sweep on mass and the number of runs is reduced further, to: 

N total N LdN return— date (13) 

The Jacobian sparsity pattern for the MSSR problem with variable pickup mass, spiral descent /ascent 
modeling, and no planetary flybys is shown in Figure 4. The entries common to all three versions of the 
problem are shown in blue, and the entries common to the fixed-mass and variable-mass problems with spiral 
descent /ascent modeling are shown in purple. The new entries unique to the variable mass problem are shown 
in red. There are ten new entries. The first seven are a vertical line in the thick blue bar representing the 
dependencies of the return-journey match point constraints. These represent the dependency of the match 
point constraints on the variable pickup mass. The final three entries are immediately below the thick blue 
bar and represent the dependence of the Earth return date and total flight time on the pickup mass. These 
dependencies are introduced because when the pickup mass changes, the time required for the ascent spiral 
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changes and influences the return date and total flight time. In addition, the pickup mass influences the 
amount of propellant necessary for the ascent spiral and therefore the dry mass constraint. 



The variable pickup mass with spiral modeling formulation of the MS SR problem is the most complex 
of the three presented here and requires a tool that can properly model the additional dependencies but 
significantly reduces the number of parametric runs required. This is the most efficient, in terms of the 
human overhead required to solve the problem, approach to designing MS SR trajectories and is therefore 
implemented in EMTG. 

III.D. Stochastic Global Search via Monotonic Basin Hopping and SNOPT 

Like most other optimization problems, low-thrust trajectory design problems require an initial guess. Such 
initial guesses are often generated using Lambert’s method, by solving for the optimal impulsive trajectory, by 
using shape-based trajectory approximations, or by experience-driven intuition. Once an initial guess is found 
for one version of the problem, the solution to the first problem may be used as an initial guess for related 
problems. For example if one finds a solution to the fixed-mass rendezvous MS SR problem with 2000 kg 
of return mass, that solution may be used as an initial guess for the same problem with 3000 kg of return 
mass. However there are drawbacks in this approach. First, the optimal low-thrust solution may not be in 
the neighborhood of the optimal solution to the low- fidelity Lambert version of the same problem. Second, 
sometimes the solution space bifurcates - that is, sometimes as the pickup mass increases multiple families 
of solutions can be found and if one uses a previous solution with lower pickup mass as an initial guess, 
one falls into an inferior solution family. Similarly, if one is solving the variable pickup mass problem, the 
optimizer will only converge to the locally optimal solution that represents the best in the solution family 
that was used as the initial guess. 

Recent research in low-thrust trajectory optimization has led to the creation of stochastic search methods 
that do not require an initial guess. 22, 23, 24, 25, 26, 6 ’ 7 ’ 8 ’ 9 * ’ 10, n ’ 12 Such methods are ideal for the MSSR problem 
because they can be used to directly solve the low-thrust problem without the risk of introducing a bias 
by using a previous solution to a related problem or a low-fidelity approximation to the same problem. In 
addition, because each problem instance is now independent of all of the others, stochastic search methods 

9 of 22 


American Institute of Aeronautics and Astronautics 


allow many different versions of the MSSR problem to be solved simultaneously in parallel. The stochastic 
search method used in this work is monotonic basin hopping (MBH). 

MBH 27 is an algorithm for finding globally optimal solutions to problems with many local optima. MBH 
works on the principle that many real-world problems have a structure where individual local optima, or 
“basins” tend to cluster together into “funnels” where one local optimum is better than the rest. A problem 
may have several such funnels. MBH was originally developed to solve molecular conformation problems 
in computational chemistry, but has been demonstrated to be effective on various types of interplanetary 
trajectory problems. 26, 28, 29, 7 ’ 8 ’ 9 

First, an initial point x is randomly chosen. The NLP solver is run using x as the initial guess. If the 
NLP solver finds a feasible solution, then that new point x* is adopted as the new current point. If the 
NLP solver does not find a feasible solution, then a new random point is chosen. Once a feasible solution 
is found, MBH will attempt to “hop” from the feasible and locally optimal x* to a better point. This is a 
two-step process: first a small random perturbation vector is added to x*, producing a new x', and then the 
NLP solver is run. If the resulting solution is both feasible and superior to x*, then it is adopted as the new 
x* and the hopping process begins again. Otherwise, MBH attempts a new hop from the current x*. Each 
feasible solution is stored in an archive. 

MBH is run until either a specified number of iterations (trial points attempted) or a maximum CPU time 
is reached, at which point the best solution stored in the archive is returned as the solution to the outer-loop. 
The version of MBH used in this work has two parameters - the stopping criterion and the type of random 
step used to generate the perturbed points x'. In this work the random step is drawn from a bi-directional 
Pareto distribution with the Pareto parameter, <a, set to 1.5. The bi-directional Pareto distribution will 
usually generate small steps that allow MBH to exploit the local funnel around the current best solution. 
However some of the steps generated by the bi-directional Pareto distribution will be much larger, in some 
cases spanning the entire solution space. These larger steps allow MBH to explore the full problem. This 
approach is known to be robust on complex low-thrust problems. 12 

In addition, SNOPT occasionally freezes in the middle of a local optimization. This behavior disrupts 
the MBH global search, so EMTG contains a timer that ends any SNOPT run that continues for longer than 
some threshold time typically set to a few minutes. The pseudocode for MBH is listed in Algorithm 1. 

Algorithm 1 Monotonic Basin Hopping (MBH) 

while not hit stop criterion do 
generate random point x 

run NLP solver to find point x* using initial guess x 
if (x* is a feasible point) then 

X current — X 

while A',,,,! improve Max not 

improve do 

generate x' by randomly perturbing x current 

run NLP solver to find point x* using initial guess x' 

if x* is feasible and /(x*) < /(x current ) then 

X current — X 

end if 

end while 

save x* to archive 

end if 
end while 

return best x* in archive 


The MBH+ SNOPT optimization algorithm in EMTG is efficient and does not require an initial guess. 
MBH is most useful when one does not have much a priori information about the solution as is often the 
case in the MSSR problem. In particular, MBH can find the optimal solution to the variable-mass with 
spiral descent /ascent modeling MSSR problem when there is no initial guess for the pickup mass (or any 
other variable in the problem). The combination of the variable mass with spiral descent /ascent modeling 
formulation with the MBH+SNOPT optimization algorithm is therefore a very valuable tool for MSSR. 
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IV. Results 


IV. A. Overview 

Trajectories are designed to each of the three target orbit about Mars using EMTG. Both the FH and SLS 
departures are evaluated in addition to the 2000 and 3000 s I sp spacecraft configurations. While the the 2000 
second configuration increases the spacecraft dry mass and is less efficient than the nominal 3000 second 
configuration, it provides higher thrust. This additional thrust at the cost of efficiency is often beneficial 
in short time-of-flight scenarios in which the spacecraft thrusts nearly continuously. In these scenarios the 
maximum return mass can be limited by the thrust capability of the engine and not the propellant of the 
spacecraft. For trajectories generated with the higher I sp system, it can be more optimal to depart Earth 
with less than the maximum allowable propellant. 

A range of launch and departure opportunities is inspected, broadening the trade space and allowing for 
more informed mission planning. The search for optimal return-mass trajectories is segmented into one year 
launch and one year arrival periods. The one year launch periods for Falcon Heavy departure begin on the 
September 23rd of 2018, 2019, 2020, 2021, or 2022. The one-year arrival period spans the calendar year 
from January 1 to December 31, where 2024, 2025, 2026, and 2027 arrival years are investigated. For SLS 
launches the one- year launch date period begins on July 1. These maximum allowable return mass matrices 
can then be used to determine when cost and benefit of different launch and arrival date combinations. For 
certain Earth departure and arrival date periods, gravity assists can improve the maximum allowable return 
mass. Gravity assists from Earth and Venus are considered for comparison to direct trajectories. Earth 
gravity assists on ARRM Option B return legs have demonstrated the potential to significantly increase the 
maximum boulder return mass for some of the NEAs examined, 2 and are especially beneficial for changing 
the trajectory plane from the inclination of the asteroid to the ecliptic for Earth-Moon system insertion. For 
the outbound leg, both Earth and Venus are explored, where gravity assists from those bodies have been 
shown to improve low-thrust trajectory performance on outbound journeys to Mars. 30 

IV.B. Deimos Results 

Deimos orbits Mars at the highest altitude of the options investigated and, thus, provides the highest feasible 
return masses. As such, the body provides an upper bound on performance. For trajectory modeling purposes 
Deimos is assumed to reside in circular orbit about Mars at a radius of 23453 km. Another key spiral modeling 
assumption is that the spacecraft arrives at Mars sphere of influence with a C3 of zero, and thus the spiral 
can begin at any Mars inclination. The density of Deimos is 1.48 g/cra 3 , 31 and this value can be used to 
estimate the size of the boulder size that can be returned from the pickup mass. The stay time at Deimos 
for proximity operations is constrained to be greater than 200 days, but less than 600 days (also used for 
Phobos and LMO trajectories). 

As noted, different interplanetary options including gravity assists to and from Deimos are explored, 
namely: EDE (direct), EEDE, EDEE, EVDE, EEDEE, EVDEE, EVEDE, and EVEDEE sequences, where E, 
D, and V indicate an encounter with Earth (including departure and arrival), Deimos, or Venus respectively. 
The single Earth gravity assist (outbound or inbound) options are generally more beneficial than trajectories 
with a Venus gravity assist or a trajectories with a gravity assist on both legs. Trajectory results are listed 
in Table 2 for a one year launch periods beginning no earlier than (NET) September 23, 2019 and a one year 
arrival period beginning no later than (NLT) January 1, 2027. For these particular launch and arrival periods, 
the direct outbound and return leg trajectory outperforms any of the gravity assist scenarios. The return 
mass from the EEDE sequence is, however, only 200 kg lower than the direct scenario. The scenarios using 
a Venus gravity assist do not perform particularly well given the long outbound time-of- flights. High return 
mass solutions typically benefit from allowing as much flight time as possible for more optimal thrust /coast 
periods and phasing. Note that the EVDEE and EVEDEE scenarios are either not feasible or return less than 
460 kg of asteroid mass. The trajectories are identified as infeasible in EMTG; however, 460 kg of hydrazine 
and xenon is reserved as part of the neutral mass for proximity operations at Deimos. This reserved mass 
does not need to be returned, and the maximum return mass identified using EMTG can be augmented by 
460 kg. Thus, the scenarios may be feasible, but are not competitive in maximum return mass. 

Understanding the performance of trajectories over different combinations of departure and arrival dates 
is of critical importance for the planning of Mars moon or low-Mars orbit sample return mission. The 
maximum Deimos return mass with different gravity assist scenarios and I sp configurations over the different 
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Table 2: Trajectory characteristics and performance for direct and gravity assist scenarios departing Earth 
NET 9/23/19 (one year period) and returning NLT December 31, 2027 (one year period) 


Scenario 

Earth System 
Departure 
Date 

Demos 

ArrhalDate 

Stay 

Time 

(days) 

Deimos 

Departure 

Date 

Earth Skstem 
ArrhalDate 

TOF 

(years) 

Launch 

Mass 

(kg) 

Asteroid 
Arrival 
Mass (kg) 

Earth 

Return 

Mass (ks) 

Boulder 

Mass 

(tons) 

EDE 

7-Jait-2020 

'-Oct-2022 

200.0 

4192023 

3 1-Qc t-202 7 

7. 32 

12242 

9610 

15315 

9.937 

EEDE 

j-Jaur-2020 

17-Aus-2022 

269.4 

13-Nfey-2023 

2S-Oct-2027 

7.32 

12024 

9335 

15093 

9.737 

EDEE 

29-Sei>2019 

27-Nir-2022 

200.0 

I3-Oct-2022 

3 -Feb-2027 

7.36 

:0919 

3231 

13281 

3.091 

EYDE 

20- Nov-2019 

9-JurP2024 

200.0 

26-Dec-2024 

S-Oct-2027 

7.33 

12310 

3295 

11334 

6462 

EEDEE 

1B-Oct-2019 

2 3 -Feb-2022 

200.0 

16-Sep-2022 

29-Jafr202.7 

7.23 

10190 

3214 

13329 

3.710 

E\DEE 

Not feasible or return mass. <460 ks 

EVEDE 

23- Sep-2019 |25-Nhy-2024| 

1 200.0 | 

1 11 -Dec-2.024 25 -Nov- 2027 

1 3.17 

1 11293 | 

1 7397 

1 11046 

1 5.326 

EYEDEE 

Not fea able or return nia ss < 460 kt 


Earth system departure and arrival periods under consideration are arrayed in Tables 3-5. A in the 
tables indicates that there are no feasible solutions or the return mass is less than 460 kg for that date 
combination. The entries are also color coded according to their relative performance, where the darker the 
green the higher the maximum return mass. Table 3 outlines the maximum return mass for direct, EDE 
scenarios for both the 2000 second and 3000 second I sp configuration. The 3000 second I sp configuration 
generally outperforms the 2000 second I sp configuration except for the 2024 arrival dates, which correspond 
to short time of flight (TOF) trajectories that benefit from the higher thrust for improved phasing. The 
EEDE and EDEE trajectories in Tables 4 and 5 follow this same trend, with the 3000 second I sp configuration 
usually enabling higher return masses except when the TOF is relatively short. 


Table 3: Maximum return mass for EDE trajectories with a FH launch and different I sp configurations 


Sequence: EDE . Launch vehicle: FH. Isp: 2000s Sequence: EDE. Laiuch vehicle: FH. Isp: 3100s 


Earth system return war 


Earth system return year 


>■ 

^ P H=- 

p f-' 

E ^ 

- | ™ 

^ — 



2024 

2025 

2026 

2027 



2024 

2025 

2026 

2027 

20 IS 

1.278 

7.124 

6.648 

7.143 

i 

2018 

0.540 

9.410 

7.413 

10.413 

2015 

1.121 

7.169 

6.652 

7.189 

g a 2 

2019 

0.551 

7.839 

6.876 

9.987 

2020 

- 

0.695 

- 

7.042 

a lu 

— :D K'i 
“D Eu 

2020 

- 

- 

- 

8.685 

2021 

- 

- 

- 

7.093 

Z 

2021 

- 

- 

- 

7,675 


Table 4: Maximum return mass for EEDE trajectories with a FH launch and different I sp configurations 


Sequence: EEDE , Launch vehicle: FTL Isp: 2000s Sequence: EEDE., Launch vehicle: FHTsp: 3 0 00s 


Earth system retun year 


Earth system retun year 


P 

I s.F 

I I % 

= % S 

p] Jf r- 
J & 



2024 

2025 

2026 

2027 

2018 

1.936 

8.176 

7.967 

■- j 
--j 

s 

2019 

1.932 

7.206 

7.122 

3.566 

2020 

- 

- 

- 

- 

2021 

- 

- 

- 

7.980 


& 9 > 

K ^ >- 

g S 

JS -H 

B a ^ 


& 



2024 

2025 

2026 

2027 

2018 

0.714 

3.321 

6,990 

9.557 

2019 

- 

3.315 

6.634 

9.600 

2020 

- 

- 

- 

- 

2021 

- 

- 

- 

3.479 


A mixture of sequences and I sp configurations comprise the best performing cases across the different 
departure/arrival date combinations. The best Deimos-FH trajectories in terms of return mass for each of 
the date combinations are compiled in Table 6 and the corresponding trajectory details are listed in Table 
7. All trajectories depart with the maximum possible C3 of 2 km 2 /s 2 . Short TOF trajectories (<5.3 years) 
benefit from the higher thrust of the 2000 second I sp configuration. Note that the launch mass is at its 
upper bound of 13153 kg for the 2000 second I sp configuration trajectories to compensate for the higher 
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Table 5: Maximum return mass for EDEE trajectories with a FH launch and different I sp configurations 


Sequence: E DEE lunch vehicle: FH, Isp: 2000s Sequence: E DEE . Launch vehicle: FH. Isp: 3000s 


3 S ^ 

-j -g *3 
t] a ch 

^ J" b 
" & 


Earth system return year 
(XLT 12./3 1 /%%%%) 


Earth system return year 
(NTT 12.31 Axes) 



2024 

2025 

2026 

2027 

r- b ^ 


2024 

2025 

2026 

2027 

2018 

- 

4.625 

7.949 

7.954 

i tt 

* o £; 
f ' g 

- -g 

2018 

- 

4435 

9_S26 

9.860 

2015 

- 

- 

7,848 

7.934 

2015 

- 

4.405 

8.061 

8. 101 

2020 

- 

- 

- 

3. 838 

T“ n Ch 

£ f 5 

2020 

- 

- 

- 

3.260 

2021 

- 

- 

- 

4.312 

e. 

2021 

- 

- 

- 

2_5 86 


propellant burn rate. On the other hand, optimal trajectories with a 3000 second I sp configuration typically 
depart with less than the maximum allowable propellant mass to avoid returning to the Earth system with 
any unused xenon. Notably, an Earth gravity assist on the outbound leg is beneficial for several of these 
short TOF cases, improving the return mass over the direct scenarios in some departure/arrival date pairs 
(2018/2024, 2019/2024, 2019/2025, and 2021/2027). With a longer TOF, an Earth gravity assist on the 
return leg can reduce the propellant required to arrive at the Earth system with a C3 of 2 km 2 / s 2 . An 
EDEE sequence with a 3000 second I sp is the highest performing scenario for the 2018/2026 and 2019/2026 
departure/arrival date combination. 


Table 6: Best return mass as a function of departure and arrival period across all Deimos trajectories 

examined with a FH launch 
Earth system return war 


i if 
1 1 s 
! 3 a 
- # U 

z 



2024 

2025 

2026 

2027 

201S 

1.936 

9.410 

9.826 

10.413 

2015 

1.932 

8315 

8.061 

9.987 

2020 

- 

0.695 

- 

8.685 

2021 

- 

- 

- 

8.479 


Table 7: Trajectory characteristics of best return mass cases across all Deimos trajectories examined with a 

FH launch 


E arth System 
Departure Date 

E arth System 
Arrival Date 

Sequence 

Isp 

(sec) 

Deimos 
Arrival Date 

Stay 

time 

(days) 

TOF 

(days) 

L aimch 

Mass 

(Lsj 

Deimos 
Arrival 
M ass 

(to 

Earth 

Return 

Mass 

Q& 

Max 

Boulder 

Mass 

(tons) 

20-_%e-2019 

3 1 -Dec-2024 

EEDE 

2000 

16- Dec-2021 

200 

5.37 

13153 

9S92 

7475 

1936 

1 6- Dec- 201 S 

'9- Sep-202 5 

EDE 

3000 

ll-Jkn,2022 

200 

6 76 

12220 

S7SI 

14734 

9410 

18-Dec-201S 

7 -Oct- 2026 

EDEE 

3000 

30- Dec-2021 

200 

7. SO 

12067 

S664 

15136 

9.S26 

23- Sep-2019 

6-Dec- 2027 

EDE 

3000 

4- Jut 2022 

200 

S.20 

13096 

10109 

15S27 

10413 

26- Sep-2019 

3 0-Dec- 2024 

EEDE 

2000 

21-Dec- 2021 

200 

5.26 

13153 

9S73 

7472 

1932 

2 S- Sep-2019 

20- Sep-202 5 

EEDE 

3000 

14-\fei-2022 

200 

5.9S 

10244 

S3 16 

13439 

S.315 

2 S- Sep-2019 

9-Oct-2026 

EDEE 

3000 

2S-Mu-2022 

200 

7.03 

10967 

S250 

13255 

S.061 

7-Jau,2020 

3I-Oct-2027 

EDE 

3000 

I-Oct-2022 

200 

7.S2 

12242 

9610 

15315 

9.9S7 

4- Dec-2020 

l-Jun,2Q25 

EDE 

2000 

2- Jut 2023 

200 

449 

13153 

7956 

622 S 

0695 

IS- Jan, 2021 

25-'Oct-2027 

EDE 

3000 

20-Feb-2024 

200 

6 76 

12299 

S737 

14014 

S6S5 

2S-Dec-2021 

19-Oct-2027 

EEDE 

3000 

24-yfer-2024 

200 

5.81 

10906 

S513 

13669 

S.479 


The best 2019/2025 trajectory is an EEDE sequence with a 3000 second I sp configuration. That trajectory 
is plotted in Figure 5, illustrating an energy-raising Earth gravity assist within one year of departure after 
a period of leverage thrusting. After the flyby, the spacecraft then thrusts continuously (90% duty cycle) 
to reach Mars, arriving at the Mars SOI 15.5 months later with a C3 of zero. The spiral down to Deimos 
orbit requires approximately 1.08 km/s 2 of Av and five months of flight time with starting mass of 8626 kg. 
The stay time for proximity operations at Deimos is 200 days, during which the spacecraft can retrieve up 
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to 8.3 tons from Deimos before continuously thrusting to spiral back out to Mars SOI and back to the Earth 
system. 


Ie8 


x (km) 


-3 

- 4 — 


-2 


-1 

- 4 - 



■“ -2 


- 2 , 


Figure 5: Ecliptic projection of highest Deimos return mass trajectory with a 2019 departure and 2025 
arrival for a FH launch (sequence: EEDE, I sp : 3000 seconds) 


Optimal return-mass trajectories for Deimos are also examined for SLS launches. An Earth gravity assist 
on the return leg (EDEE) is considered in addition to direct sequences (EDE). Given the high capability 
of the launch vehicle, a gravity assist on the outbound leg are not examined. Direct and gravity assist 
return masses from Deimos as a function of one- year Earth launch and arrival date period are outlined in 
Table 8. The one year launch period begins on July 1 for the SLS cases, and the 2019 launch period cases 
are not shown as the optimal trajectories are effectively the same trajectory as the 2020 trajectories for 
corresponding arrival dates. All SLS trajectories are constrained so that the first 30 days of the trajectory 
is a coasting period, allowing for check-out of the SEP system. Only return masses for the 2000 second I sp 
configuration are shown as they are higher than the return masses for the 3000 second configuration in all 
launch and arrival date combinations considered. With the additional launch capability, the SLS is able to 
deliver more propellant to a higher C3 and the higher thrust of the 2000 second I sp configuration can be 
leveraged. 

The characteristics of the best performing SLS trajectories are listed in Table 9. In only the longest 
TOF launch/arrival date pair (2020/2027) is it optimal to launch the maximum spacecraft mass of 16072 
kg. This trajectory is the highest performing out of all launch/arrival date pairs with a return mass of over 
13 tons, and is plotted in Figure 6. The other scenarios are constrained by the shorter TOF and launching 
less propellant is optimal. Unlike the FH launch trajectories, gravity assists are not as beneficial for SLS 
launches, where the EDEE gravity assist trajectory out-per forms the direct trajectories in only one date 
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combination (2020/2026). 

Table 8: Maximum Deimos return mass as a function of departure and arrival period for EDE and EDEE, 

2000 second I sp trajectories with a SLS launch 


Sequence.' EDE, Launch vehicle: SLS, Isp: 2000s Se£|ljeiKe . E D EE, Lamch vehicle: SLS, Isp: 2000s 


K *" 1 

9 > 

.-j - 

Tj 

j= r-" 




Earth B^tem return year 
(NLT 1231/YYYY) 



2024 

2025 

2026 

2027 

2020 

2.097 

10.661 

9.049 

13.454 

2021 

- 

- 

- 

11.607 

2022 

- 


- 

11.740 


E ar±L system return war 


> 
j E 

j= r~- 



2024 

2025 

2026 

2027 

2020 

- 

6.165 

11.157 

11.272 

2021 

- 

- 

- 

4.869 

2022 

- 

- 

- 

5.091 


Table 9: Trajectory characteristics of best Deimos return mass cases across all departure and arrival date 

combinations examined with a SLS launch 


Earth System 
Departure 
Date 

E arth System 
Arrival Date 

Sequence 

i— i 

a 

(km A 2/s A 2) 

Deimos 
Arrival Date 

Stay 

lime 

(days) 

TOF 

(days) 

Launch 

Mass 

fkE) 

Deimos 

Arrisral 

Mass 

(ke) 

Earth 

Return 

Mass 

CkE) 

Max 

Boulder 

Mass 

(tons) 

1 -Jut 2020 

3I-Dec-2024 

EEE 

2000 

S.49 

17 -Dec-2021 

200 

450 

122 IS 

9S92 

7531 

2.097 

4-Jut 2020 

20^ Sep-2025 

EDE 

2000 

9.30 

20-Apr-2022 

200 

5 21 

14217 

S7S1 

16323 

10.661 

1 -Jut 2020 

30-Dec-2026 

EDEE 

2000 

9.00 

4- Apr- 2022 

200 

6.50 

14004 

S664 

16799 

11 157 

2S-Jut2020 

24- Oct-202 7 

EEE 

2000 

14.5S 

3 -Sep- 2023 

200 

7.24 

16072 

10109 

19310 

13 454 

30-Juit-2022 

6- Oct-202 7 

EEE 

2000 

12 19 

10- Apr- 2024 

200 

5.27 

14917 

9S73 

17342 

11607 

21-Jut 2022 

7-Gct-2027 

EEE 

2000 

1015 

9- Apr- 2024 

200 

5 21 

14S33 

11SS6 

17467 

11.740 


IV. C. Phobos Results 

A similar search for high Phobos return mass scenarios across different departure and arrival date combi- 
nations is also conducted. The key difference in the design is the additional time and propellant required 
to spiral down to Phobos orbit radius of 9234 km. The uniform density of Phobos is estimated to be 1.87 
g/cra 3 , 31 which is slightly higher than Deimos. 

As with Deimos, both direct and gravity assist trajectories are designed for 2000 second and 3000 second 
I sp configurations. EPE, EEPE, and EPEE scenarios are examined for FH launches, and the return mass 
from the best performing trajectories for each of the one-year departure/arrival periods is outlined in Table 
10. The trajectory characteristics associated with these optimal return masses are listed in Table 11. Notably, 
the sequences are the same as optimal Deimos trajectories for corresponding departure/arrival date period 
combinations except for the 2020/2027 pair. The optimal I sp configuration is, however, not the same. 
Additionally, there are not feasible trajectories that return more than 460 kg for 2018/2024, 2019/2024, and 
2020/2025 pairs because of the additional time and propellant required for spiraling (identified in Table 10 
by Substantial return masses of more than 9 tons are feasible, however, with a departure in 2018 or 

2019 and an Earth return in 2027. These two entries are effectively the same trajectory with very similar 
associated dates and masses. The optimal trajectory associated 2019/2027 entry is plotted if Figure 7. In 
this trajectory the spacecraft thrusts continuously, with the spiral down to the orbit of Phobos requiring a 
Av of 1.88 km/s and nine months of flight time. 

Trajectories to Phobos are also identified for SLS launches. The results across the launch and arrival date 
periods considered are shown in Table 12. The same constraints and assumptions as for Deimos trajectories 
are applied. Only direct sequences are explored as it is less likely that a gravity assist on the return leg will 
offer improvement given the additional time required for spiraling as compared to Deimos. The 2000 second 
configuration enables higher return masses in all date combinations considered, and the associated trajectory 
characteristics are outlined in Table 13. As with the best Deimos SLS trajectories, the maximum amount of 
propellant is launched in the one scenario (2020/2027) associated with the longest TOF. The SLS launch 
vehicle allows for later departure dates and shorter TOFs as compared to the FH launch trajectories. 
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Figure 6: 


Table 10: 


leS 


x (km) 



Ecliptic projection of highest Deimos return mass trajectory with a 2020 departure and 2027 
arrival for a SLS launch (sequence: EDE, I sp : 2000 seconds) 


Best return mass as a function of departure and arrival period across all Phobos trajectories 

examined with a FH launch 


Earth system return year 
(XLT 12.-3 l.'YYYY) 



2024 

2025 

2026 

2027 

201S 

- 

5.249 

5.331 

9.13S 

201 9 

- 

3.668 

4.55 1 

9.164 

2020 

- 

- 

- 

5.407 

2021 

- 

- 

- 

4.840 
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Table 11: Trajectory characteristics of best return mass cases across all Phobos departure and arrival date 

combinations examined with a FH launch 


E arth System 
Departure Date 

E arth System 
Arrival Date 

Sequence 

Isp 

(sec) 

Phobos 

ArrivaLDate 

Stay 

time 

(days) 

EOF 

(days) 

Launch 

Mass 

(kg) 

Phobos 

Anris-al 

Mass 

(kg) 

Earth 
Return 
M ass 
(kg) 

M as 
Boulder 
Mass 
(tons) 

30-Nov-201S 

15- Sep-2 025 

EPE 

3000 

6- Apr-20 22 

200 

6. SO 

11S30 

S223 

10519 

5.249 

25-NOY-201S 

3- Sep-2026 

EPEE 

3000 

31-Mtr-2022 

200 

7.77 

11S09 

S1S3 

10599 

5.331 

22- Sep-2 019 

2 6- Oct-2 027 

EPE 

3000 

14- Oct-2 022 

200 

S.09 

12695 

9520 

1450 S 

9.137 

l-Oct-2019 

3I-Dec-2025 

EEPE 

3000 

6-Jun-2Q22 

200 

6.25 

9972 

7S73 

S74S 

3.66S 

2 S- Sep-2 019 

1 7-Nov-2026 

EPEE 

2000 

23-Ajie-2022 

200 

7.14 

13153 

959 S 

10103 

4.551 

27-Oct-2019 

2 6- Oct-2 02 7 

EPE 

3000 

1 5-Oct-2 022 

200 

S.OO 

12546 

9506 

14519 

9164 

15- Apr-2 021 

1- Sep-2027 

EEPE 

2000 

IS- Jut 20 24 

20 S 

6.3 S 

13153 

LOO 99 

10964 

5.407 

23- Sep-2 021 

:- Sep-202 7 

EEPE 

2000 

26-Jun-2024 

231 

5 94 

13153 

9SS0 

10394 

4.S40 


le8 


x (km) 


-3 


-2 


-1 



-1 


Figure 7: Ecliptic projection of highest Phobos return mass trajectory (2019 departure and 2027 arrival) 

for a FH launch (sequence: EPE, I sp : 3000 seconds) 
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Table 12: Maximum Phobos return mass as a function of departure and arrival period for direct 2000 and 

3000 second I sp trajectories with a SLS launch 


Sequence: EFE. Launch vehicle: SLS, Isp: 2000s 
Earth system return ye ar 


Sequence: EFE, Launch vehicle: SLS. Isp: 3000s 
Earth system return ye ar 


(NLT 12/3 1/YYYY) (XLI 12/31/YYYY) 



2024 

2025 

2026 

2027 



2024 

2025 

2026 

2027 

2020 

- 

6.057 

5.074 

10 605 

it 

2020 

- 

4.599 

3 .503 

9.252 

2021 

- 

- 

- 

6.625 

•Jf r_ 

2021 

- 

- 

- 

- 

2022 

- 

- 

- 

6.635 

—■ 1ZJ 
“ & 

2022 

- 

- 

- 

4.7 8 5 


Table 13: Characteristics of maximum Phobos return mass trajectories with a 2000 second I sp 

configuration and a SLS launch 


E arch System 
Departure Date 

Earth System 
ArrivalDate 

Sequence 

Isp 

(sec) 

o 

(km A 2/s A 2) 

Phobos 

ArrivalDate 

Stay 

tine 

(days) 

TOF 

(days) 

Launch 

Mass 

(tg) 

Phobos 

Arrival 

Mass 

m 

Earth 

Return 

Mass 

m 

Mat 

BoiMder 

Mass 

(tons) 

1-M2020 

19-Aie-2025 

EPE 

2000 

8.73 

20-Juht-2022 

200 

5. 14 

13575 

10417 

11639 

6057 

l-M-2020 

1-J.aiv2026 

EPE 

2000 

S.44 

25-Jfiv2(E2 

200 

5 JO 

131S3 

10142 

10614 

5.074 

22-Jut2020 

2^Oct-2027 

EPE 

2000 

13.87 

23 -Nov- 2023 

200 

7.27 

16072 

123 S9 

16451 

10605 

30-Ji^2022 

30-^-2027 

EPE 

2000 

1403 

23-Juor2024 

200 

5.17 

14152 

10731 

12287 

6.625 

S-Jut2022 

30-Ale -202 7 

EPE 

2000 

13 47 

23-Junr2Q24 

200 

5.15 

14094 

10729 

12292 

6.635 


IV. D. LMO Results 

Trajectories are additionally designed to evaluate the ability of the ARRV to spiral down to a 400 km altitude 
Mars orbit to retrieve a Mars sample. The same driving assumptions and constraints as in the Deimos and 
Phobos cases are made, where the 460 kg of propellant reserved for landing on one the moons is instead used 
for rendezvous with the sample capsule. 

For FH launch scenarios EME, EMEE, and EEME (M indicates Mars for LMO) sequences and 2000 and 
3000 second I sp configurations are surveyed. The best return masses from LMO for FH launches as a function 
of Earth system departure and arrival date period are arrayed in Table 14, and the associated trajectory 
characteristics are in Table 15. As with Deimos and Phobos, an assortment of different sequences and I sp 
configurations are represented for FH launches. Gravity assists are again advantageous for improving the 
maximum return mass, in which five of the eight best trajectories exploit an Earth flyby on the outbound or 
return leg. Similar trends are also observed on the best I sp configuration for a given departure/arrival period 
entry, where a 2000 second I sp configuration is better for the shorter TOF cases and a 3000 second I sp config- 
uration is better for the longer TOF scenarios. An early launch and a relatively long TOF enable maximum 
return masses greater than six tons. If an additional 1700 kg payload is carried to LMO with the ARRV, 4.8 
tons can still be returned from LMO to an Earth-return C3 of 2.0 km 2 / s 2 for the 2018/2027 and 2019/2027 
date combinations. This additional payload could represent a Mars surface descent /retrieval/ascent system, 
enabling a Mars sample return with a single FH launch. 


Table 14: 


Best return mass as a function of departure and arrival period across all LMO trajectories 

examined with a FH launch 


Earth ostein return year 
CN'LT 12/31/YYYY) 



2024 

2025 

2026 

2027 

20 IS 

- 

1.280 

1.323 

6.225 

2019 

- 

1.101 

0.918 

6.224 

2020 

... 

- 

- 

0.588 

2021 

- 

- 

- 

0.889 
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Table 15: Trajectory characteristics of best return mass cases across all LMO departure and arrival date 

combinations examined with a FH launch 


E arth System 
Departure D ate 

E arth System 
Arrival Date 

Sequence 

Isp 

(sec) 

LMO Arrival 
Date 

Stay 

time 

(days) 

TOF 

(days) 

Launch 

Mass 

(kg) 

LMO 

Arrhal 

Mass 

(kg) 

Earth 

Return 

Mass 

(kg 

Max 

Boulder 

Mass 

(tons) 

16-Xlav2019 

22-Ju>2 025 

EEME 

2000 

3- Sep-2022 

200 

6.10 

13153 

9116 

6316 

1.230 

14-Nov^20l3 

2 5 -Oct- 2026 

EMEE 

3000 

i-Jii-2022 

200 

7.95 

11260 

7405 

6523 

1.323 

23- Sep-201 9 

6 -Gct-2027 

EEME 

3000 

2-Jii-2023 

200 

S.04 

12277 

9192 

11543 

6.225 

23- Sep-201 9 

23-Ju>2025 

EEME 

2000 

17- Sep- 20 22 

200 

5.75 

13153 

9022 

6636 

1.101 

2 7- Sep-201 9 

3 -Oct-2 026 

EMEE 

2000 

16- Oct- 20 22 

200 

7.03 

13153 

S639 

6452 

0.91 S 

lS-Xanr2020 

6 -Oct-2027 

EME 

3000 

4-1x1-202 3 

200 

7.72 

126S4 

9244 

11533 

6.224 

2 3- Sep-2021 

7-Jl42027 

EME 

2000 

14-N ov-2024 

200 

5.79 

13153 

3309 

6120 

0.533 

3-JiL2022 

S-Jl42027 

EME 

2000 

22-Oct-2024 

200 

5.01 

13153 

S93 S 

6423 

0.339 


Returns from LMO are viable for TOFs as short as five years with a FH launch. The best 2021/2027 
trajectory has a 5.01 year TOF and is plotted in Figure 8. Interestingly, the outbound and return leg legs 
are nearly symmetric in shape with Mars line of apsides. This geometry enables the spiraling period down 
to LMO and back to Mars SOI to occur around Mars perihelion instead of aphelion, allowing for higher 
power to the SEP system. The spiral down to LMO requires a Av of 3.09 km/s and 10.5 months. After the 
sample is retrieved the spiral back up to SOI takes 11.5 months. Recall that the arrays degrade throughout 
the trajectory decreasing the power available. 

SLS launch scenarios to LMO are designed as well. As with Phobos, only direct sequences are examined 
for both I sp configurations. The maximum return mass as a function of launch and Earth system arrival 
date period are illustrated in Table 16. The trajectory characteristics of the best performing cases for 
the launch/arrival periods are listed in Table 17. Unlike the Phobos SLS cases there is a 3000 second I sp 
configuration scenario that enables a higher return mass than the 2000 second I sp configuration. Returns are 
feasible for the same launch/arrival date periods as Phobos SLS cases, but are significantly lower. Despite 
the long spiral times, over six tons can be returned from LMO for a seven year TOF. 


Table 16: Maximum LMO return mass as a function of departure and arrival period for EME, 2000 and 

3000 second I sp trajectories with a SLS launch 


Sequence: EME. Launch rehide: SLS. Isp: 20003 
Earth system return ar 


o 



(NLT 1 2/3 1 : Y Y YY) 



2024 

2025 

2026 

2027 

2020 

- 

1.426 

0.760 

5.620 

2021 

- 

- 

- 

1 313 

2022 

- 

- 

- 

1 313 


Sequence: EME .Launch rehide: SLS.Isp: 3000s 


Earth system r eturi ye ar 
(XII 12/31 YYYYj 



2024 

2025 

2026 

202 7 

2020 

- 

0.562 

- 

6.406 

2021 

- 

- 

- 

- 

2022 

- 

- 

- 

- 





Table 17: Trajectory characteristics of best LMOreturn mass cases across all departure and arrival date 

combinations examined with a SLS launch 


Earth System 
Departure Date 

Earth System 
Arrival Date 

Sequence 

Isp 

(sec) 

a 

(hm A 2/s A 2) 

LMO Arrival 
Date 

Stay 

dme 

(days) 

TCF 

(days) 

Launch 

Mass 

m 

LMO 

Arrival 

Mass 

(ks) 

Earth 

Return 

Mass 

(ke) 

Max 

BeuHer 

Mass 

(tens) 

l-Jii-2020 

7-Juflr2025 

EXE 

2000 

3.47 

3- Sep- 2022 

200 

493 

12553 

9100 

6903 

1.426 

l-Jii-2020 

'.-Juflr-2026 

EXE 

2000 

3.43 

23-Aus-2022 

200 

.5.92 

12323 

3934 

6211 

0.760 

21-Ale -2020 

7-Oct-2027 

EXE 

3000 

12.22 

3-JuL2023 

200 

7.12 

11543 

9151 

11653 

6.406 

19- Jut 2022. 

9-Jul2027 

EXE 

2000 

20.23 

26-5 ep- 2024 

200 

5.06 

12901 

9130 

6323 

1313 

l-Jii-2022 

IO-Jidr-2027 

EXE 

2000 

19.37 

26-5 ep- 2024 

200 

5.03 

12355 

9126 

6320 

1313 
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Figure 8: Ecliptic projection of shortest time- of- flight trajectory from best LMO return mass trajectories 
(2021 departure and 2027 arrival) for a FH launch (sequence: EME, I sp : 2000 seconds) 
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V. Conclusion 


V.A. Summary 

In this work, a design technique was developed for the MSSR mission which allows the mission design tool to 
model operations at Mars and directly optimize the return mass without needing a large number of parametric 
runs. This technique is highly efficient and significantly reduces the human work-load necessary to design such 
missions. The algorithms presented in this work have been implemented in Goddard’s low-thrust trajectory 
design tool EMTG, which is available open-source from http://sourceforge.net/projects/emtg/. 

The technique developed here was then used to design several different variants of MSSR, including 
missions to Deimos, Phobos, and LMO using two different SEP systems and two different launch vehicles 
- FH and SLS. Many viable MSSR trajectories were found with launch dates between 2018 and 2021 and 
return to cislunar space between 2025 and 2027. These results show that it is possible to return as much as 
13 tons of material from Deimos, 11 tons from Phobos, and 6 tons from LMO. In several scenarios an Earth 
gravity assist on the outbound or return leg can increase the return mass from trajectories following a direct 
sequence for the same launch and arrival period. The ARRV is shown to be a highly flexible vehicle, capable 
of accomplishing the MSSR mission as well as the ARRM Option B mission for which it was originally 
designed. 

V.B. Future Work 

As the results of this study show, there is a broad space of possible trajectory options for MSSR and 
also many possible combinations of the system parameters such as propulsion and launch vehicle options. 
The problem space might be further expanded by investigating several different solar array sizes. The 
full-problem trajectory optimization technique introduced in this work significantly reduces the work load 
required to evaluate a given problem but there are opportunities for further automation that will allow the 
problem space to be more thoroughly explored in a short period of time. 

The full MSSR problem involves trading over launch year, arrival year, launch vehicle, propulsion system, 
solar array size, possible planetary flybys, and return mass to produce a non-dominated hyper-surface among 
these parameters. Future work by these authors will investigate fully autonomous methods of finding that 
surface. 
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